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1.0  INTRODUCTION 


This  r sport  addresses  several  aspeets  of  the  remote  sounding 
problem  with  respeet  to  the  use  of  SSH-1  (Special  Sensor  H)  infrared 
satellite  data.  The  SSH-1  instrnnent  was  oarried  on  four  Defense 
Meteorological  Satellite  Program  (DKSP)  satellites  and  vas  operational 
from  1976  to  1980.  Daring  that  time  SSH-1  data  vas  intended  to  be  need 
operationally  in  support  of  veather  analysis  and  forecasting  at  Air 
Force  Global  Veather  Central  (AF8VC).  However,  several  problems  with 
both  the  data  and  retrieval  systema  have  eanaed  serious  deficiencies  in 
the  use  of  the  SSH-1  soundings. 

Vith  the  future  launeh  of  SSH-2  (modified  SSH-1)  instruments,  some 
filter  changes  vill  be  instituted.  However,  the  instrument  will  still 
be  similar  to  the  SSH-1  design.  The  lack  of  suocesa  in  the  use  of  SSH-1 
data  must,  therefore,  be  addressed. 

The  oontraet  under  which  this  work  was  completed  eovers  five  (5) 
main  areas  where  help  is  needed  and  can  be  provided.  The  first  two  of 
these  objectives  are  addressed  in  this  report.  They  are  briefly: 

1.  The  assessment  of  current  statistical  retrieval  sohemes  as 

used  by  AFGVC  and  applied  to  synoptio  scale  data. 

Assistance  vith  implementation  of  new  techniques  and 
suggestions  as  to  development  or  modification  of  new 
techniques  to  be  applied  as  needed. 

2.  The  assessment  of  the  cloud-contamination  problem,  especially 

in  statistical  retrieval  methods.  This  includes  comparisons 
vith  current  cloud  elimination  techniques  as  applied  to 
infrared  data.  Suoh  technique a  are  of  vital  importance 
unless  enough  sounding  resolution  is  available  to  obtain 
seme  clear  column  radiances. 


Several  topics  related  to  these  two  objectives  will  be  disoussed  in 
the  following  sections. 

Seetioa  2  deals  with  tie  SSH-1  radianoes.  Serious  consideration 
was  given  to  the  satellite  measurements*  so  as  not  to  overlook  any  basic 
deficiencies  which  night  exist.  The  SSH  radiances  were  subjected  to 
statistical  tests  whioh  uncovered  sons  previously  known  problems.  such 
as  the  obstruction  in  the  field-of-view,  A  statistical  structure 
function  analysis  also  allowed  an  estimation  of  the  radiance  noise 
levels.  The  noise  levels  for  many  channels  compared  favorably  w  the 
predicted  noise  levels.  This  signifies  that  the  SSH-1  aeasuraae  c»  :e 
of  the  expected  quality,  and  any  deficiency  in  retrieval  accuracy  is  not 
due  to  the  radiances  themselves,  although  the  selection  of  the  ohaanel 
filter  responses  may  not  be  optimal . 

Seotiom  3  oovexs  an  important  aspeet  of  any  satellite  data  which 
are  obtained  from  scanning  radiometers.  The  changing  viewing  angles 
cause  an  extra  variation  in  the  measurements  due  to  changing  atmospheric 
path  lengths,  especially  for  those  channels  vhieh  are  least  transparent 
to  the  atmosphere.  One  problem  with  the  previously  used  senith  angle 
corrections  was  uncovered  and  a  solution  is  proposed. 

In  Section  4  the  cloud  problem  is  addressed.  Two  methods  are 
presented  to  deal  with  clouds  by  simple  means.  One  method  uses  simple 
xadianoe  thresholds  to  deteot  cloudy  values.  The  other  method  uses  the 
different  speotral  qualities  of  the  two  window  channels  whioh  will  be 
available  on  SSH- 2  to  determine  both  cloud  amount  and  temperature  (or 
height).  These  techniques  are  needed  mainly  to  eliminate  eloudy  values 
from  consideration  in  the  retrieval  process  because  infrared 
measurements  are  severely  limited  in  eloudy  situations. 


Section  5  assesses  a  regression  retrieval  scheme  as  applied  to 
SSH-1  for  both  temperatures  and  moisture.  Testing  vas  done  on  both 
dependent  data  sets  (those  used  to  develop  the  regression  coefficients) 
and  independent  data  sets  (other  than  those  used  to  develop  the 
regression  coefficients).  On  independent  data  the  retrieval  'errors' 
are  greater  than  on  dependent  data,  as  expeoted.  However,  the  results 
show  that  sufficient  information  content  does  exist  in  the  SSH-1 
radiances  to  determine  temperatures  at  moat  levels,  with  large  poaitive 
explained  variances.  In  other  words,  the  'error*  variance  ia  much  less 
than  the  total  variance.  The  moisture  results  also  show  positive 
explained  variances  at  the  lowest  atmospheric  levels,  indicating  that 
moisture  information  does  reside  in  the  SSH-1  radiances  (although  the 
moisture  information  content  is  less  complete  than  the  temperature 
information).  This  is  contrary  to  unpublished  findings  of  poor  or  non¬ 
existent  moisture  retrieval  capabilities  at  AFG1C.  Their  poor  results 
may  be  linked  to  the  statistical  data  base  usee  for  regression 
coefficient  development  or  to  errors  in  the  RAOBs  or  radianees. 

Finally,  Section  6  deals  with  transmittance  software  whioh  was 
developed  for  SSH-1.  The  main  reason  for  this  software  development  was 
to  be  able  to  apply  iterative  or  minimum-information  retrieval  sohemes 
to  SSH-1  radiances.  (Such  retrieval  techniques  do  not  require  paired 
satellite-RAQB  data  as  a  statistical  base.)  The  software  combines 
techniques  from  various  sources  for  the  H^O,  OO^,  and  window  channels 
and  it  remains  to  be  thoroughly  tested.  Development  was  discontinued  in 


lien  of  siapler  net hods  for  retrieval*  sinee  both  the  iterative 
retrievala  and  transaittanoe  calculations  require  extensive  oosiputer 
resources. 


2.0  STATISTICAL  RADIANCE  ASSESSMENT 


The  iiumtat  of  SSH-1  radisnoe  quality  it  an  approach  to  the 
problaa  of  determining  aaziaaa  temperature  aad  moisture  retrieval 
capabilities.  The  real  (aot  simulated)  SSH-1  radiaaoes  tested  vere 
calibrated  at  AFG1C  bat  aaoarreoted  for  seaith  aagle  variatioas.  Cloud 
eoatamiaatioa  vas  elimiaated  through  procedures  which  will  be  discussed 
ia  the  followiag  seotioa.  Statistical  structure  fuaetioa  aaalysis  is 
then  used  to  obtaia  aa  estimate  of  the  aoise  level  for  each  radiaace 
chaaael.  These  aoise  levels  are  thea  compared  to  those  predicted  usiag 
design  specifications  aad  measuremeats  made  before  launch  of  the 
spacecraft  (Yoder,  1975;  Klein  et  si.,  1975). 

Previous  work  using  structure  function  analysis  of  aoise  levels 
(Hillger  and  Yonder  Hear,  1979)  was  done  for  Vertical  Temperature 
Profile  Radiometer  (VTPR)  radiaaoes  from  the  NOAA-4  satellite.  Those 
results  gave  estimated  rms  radiance  aoise  levels  ranging  from  the 
predicted  radiometer  sensitivities  to  values  several  times  higher  for 
channels  which  sense  lower  in  the  troposphere.  Estimated  noise  levels 
presented  here  for  the  S8H-1  radiances  show  similar  results.  Vith  the 
atmospherio  variability  eliminated  these  estimated  noise  levels  should 
be  direotly  translatable  into  temperature  retrieval  capability  assuming 
that  no  bias  exists  between  retrieved  and  aotual  temperatures,  but  that 
step  has  not  yet  been  made. 

2.1  Noise-Level  Theory 

The  idea  of  obtaining  noise  levels  of  satellite  measurements  from 


structure  function  analysis  is  considered  equivalent  to  determining 
noise  levels  from  repeated  measurements  of  the  sane  f ield-of-view  (FOV) . 
Since  the  satellite  does  not  normally  measure  the  same  FOV  more  than 
once  per  sounding  channel  it  is  difficult  to  obtain  noise-level  data 
from  such  a  set  of  measurements.  The  satellite  instrument  is 
continually  scanning  the  surface  of  the  earth  in  adjacent  FOVs.  For  the 
SSB  radiances  these  spots  are  approximately  60  km  in  diameter,  or 
equivalently  60  km  between  spot  centers  at  nadir  along  the  seen  traek. 
The  spacing,  however,  does  inorease  with  inereasing  soan  angle.  Also, 
since  the  earth's  atmosphere  and  surface  eaoh  change  in  both  temperature 
and  emittanee  characteristics  from  FOV  to  FOV  these  additional 
variations  are  introduced  on  top  of  the  expected  variations  due  to  noise 
(including  instrumental  noise,  da ta-digi tiling  noise,  and  data 
transmission  noise). 

These  atmospheric  and  surface  variations  ean  be  redueed  in  a 
statistical  sense  by  analysing  an  extremely  large  number  of 
measurements.  The  mean  ehange  from  one  FOV  to  the  next  will  then  be 
evident.  This  mean  difference  will  be  a  combination  of  both  the  real 
horisontal  variation  between  the  two  spots  and  the  noise  in  the  two 
measurements.  Therefore,  by  extrapolating  this  mean  measurement 
difference  as  the  FOV  separation  distance  goes  to  xero,  all  horisontal 
variations  are  eliminated  and  the  only  variations  that  remain  are  due 
to  noise.  With  all  atmospheric  and  surfaoe  variation  due  to  horisontal 
gradient  eliminated,  this  process  is  equated  to  obtaining  repeated 
measurements  of  the  same  FOV. 

Oandin  (1963)  was  among  the  first  to  use  structure  function 
analysis  as  an  indirect  means  of  estimating  errors  in  meteorological 


data  seta.  Hie  spatial  structure  function  is  a  measure  of  the  mean- 
squared  difference  in  measurements  as  a  function  of  their  separation 
distance.  The  analysis  can  be  accomplished  on  a  single  synoptic 
situation  in  theory,  but  that  is  true  only  if  a  sufficiently  large 
number  of  data  points  are  available,  such  as  high-density  satellite 
measurements.  In  this  study  not  only  was  the  statistical  sample  large, 
but  satellite  radiances  from  several  passes  over  two  different  regions 
of  the  earth  vere  analyzed  to  further  eliminate  any  bias  due  to  a  one- 
look  analysis. 

Purely  random  errors  lead  to  an  exaggerated  value  for  the  structure 

2 

function  by  an  amount  equal  to  a  double  mean-squared  error  (2<r  )  in  the 

e 

measurements.  Each  of  the  measurements  in  a  pair  contributes  half  of 
this  uncertainty.  Of  course,  homogeneity  and  isotropy  with  respect  to 
the  structure  function  are  assumed.  In  other  words,  the  fields  may  not 
be  homogeneous  or  isotropio,  but  the  noise  on  the  measurements  is 
assumed  to  be.  This  is  most  likely  to  be  true  in  a  large  sample  or 
variation  of  synoptic  situations.  The  sample  used  here  is  considered  to 
be  sufficiently  large  that  the  results  are  only  slightly  dependent  on 
the  sample  under  study.  Previous  results  (Hillger  and  Vonder  Haar, 

1979)  showed  how  the  structure  function  varied  with  the  synoptic 
situation,  but  the  error  analysis  here  does  not  show  the  variation 
evident  in  the  individual  synoptic  situations  since  most  of  the 
variation  due  to  the  atmosphere  and  surface  is  extrapolated  out  of  the 
individual  structure  functions. 

2.2  An  Example  of  Structure  Analysis 

The  satellite  radiances  that  were  analyzed  were  obtained  at  several 


local  times  on  5  and  6  November  1979.  The  SSH-1  radiances  were  obtained 
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from  three  MSP  satellites  (1X3536,  1X4537,  tad  1X5539)  which  take 
Miiwownti  at  approxiaately  local  9,  7,  and  11,  respectively.  Die 
areas  eoasidered  were  a  15°  latitude  by  15°  longitude  area  froa  30°  to 
45°N  and  90°  to  105°1  centered  over  western  Oklahoas  (Figure  2.1)  and  a 
20°  latitude  by  20°  longitude  area  froa  10°  to  30°N  and  145°  to  165*1 
centered  over  Hawaii.  Back  satellite  pass  contained  a  saaple  of  up  to 
approxiaately  200  radianoe  aeasureaents  for  each  channel  considered. 
Statistical  pairing  then  provided  about  400  pairs  of  adjaoent  FOVs  and 
larger  nuabers  of  pairs  of  FOVs  separated  by  greater  distances. 

Briefly,  the  paired  values  were  categorised  into  range  gates  of  60  ka 
width,  with  divisions  at  30,  90,  150,  210  ka  etc.  Paired  FOVs  would 
therefore  fall  into  all  but  the  first  range  gate  (less  than  30  ka) . 

This  range  gate  spacing  allowed  the  paired  FOVs  to  be  categorised  into 
natural  divisions.  However,  adjaoent  diagonal  paisa  and  spots  separated 
by  several  FOVs  would  fall  into  the  appropriate  range  gates,  and  the 
divisions  would  becoae  less  olesr  for  larger  separations.  Details  of 
the  structure  equations  are  provided  in  previous  work  (Hillger  and 
Yonder  Hear,  1979) .  Running  sens  kept  track  of  the  paired  values  in 
each  range  gate  and  the  separation  distanoe,  so  the  results  represent  a 
1-dimensional  structure  (aeon-squared  difference)  along  with  a  nean 
separation  distanoe  for  the  pairs  in  eaoh  range  gate. 

One  exaaple  of  a  structure  function  for  SSH-1  channel  El  at  668 
ca  1  (C02  Q  branch)  is  shown  in  Figure  2.2.  This  graph  shows  results  of 
the  structure  analysis  on  the  Hawaii  sector  outlined  above.  A  total  of 
5  satellite  passes  went  into  the  calculations  with  over  700  radiances  in 
the  analysis.  Structure  values  are  given  at  aeon  separation  distances 
of  approxiaately  60,  120,  180,  240  ka  etc.  These  axe  halfway  between 
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Figaro  2.2  An  exsaple  of  888-1  rtdliui  structure  for  okuMl  El  (668 
ca-1)  for  thi  sab-tropic  (Hawaii)  sector.  Three  carves  have  bee* 
fitted  to  the  structure  values  by  least-squares  regressioa  in  order 
to  deteraine  the  aost  probable  intereept  at  sero  separation  die- 


the  range  gate  linits.  No  paired  FOVe  were  separated  by  less  than  30  fan 
so  the  first  oecnpied  range  gate  froa  30  to  90  ka  contained  all  adjacent 
pairs  along  the  scan  line*  Adjaoent  FOVs  between  one  scan  line  and  the 
next  are  first  included  in  the  range  gate  froa  210  to  270  ka  becanse  of 
the  220  ka  aean  separation  between  scan  lines.  This  uneven  spacing 
between  along-scan  and  cross-scan  pairs  presented  soae  analysis 
probleas*  but  these  were  overcoae  by  appropriate  selection  of  the  range 
gates.  (See  Figure  2.1  for  the  SSfl-1  scan  pattern.)  The  X's  in  Figure 
2.2  are  the  actual  (aean-squared  difference)  values.  The  structure  is 
therefore  noraally  an  increasing  function  of  distance,  but  the  decrease 
at  240  ka  is  probably  due  to  the  effect  of  cross-scan  pairs  added  at 
this  distanoe.  Structure  plots  of  soae  other  channels  did  not  show  this 
feature.  The  decrease  also  aay  be  due  to  the  fact  that  the  radiances 
were  not  corrected  for  zenith  angle  variations.  Along-scan  values  would 
contan  such  variations  and  cross-scan  values  would  not.  However,  it  is 
felt  that  this  effeot  is  saall  in  oases  where  the  natural  ataospherie 
and  surface  variations  are  large,  since  those  variations  will  in  general 
be  randoa  with  respect  to  separation  direction. 

Each  structure  value  represents  the  aean-squared  difference  between 
values  at  that  separation,.  Increasing  structure  with  distanoe, 
therefore,  represents  the  aean  gradient  in  the  field  of  analyzed 
radiances.  This  gradient  is  then  eliainated  by  extrapolation  of  the 
strucuture  to  zero  separation  distance  (finding  the  intercept  with  the 
structure  axis).  This  intercept  was  found  by  the  least-squares  fitting 
of  several  lines  to  the  actual  structure  values.  The  lines  fit  to  the 
values  in  Figure  2.2  were  a  linear,  a  symetric  quadratic,  and  an 
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exponential  with  the  fossa: 


I  -  A  +  Bx 


Y  -  A  +  Ox  +  Cx  (B  -  0) 


Y  •  A  exp(Bx) 


(2.1) 


(2.2) 


(2.3) 


The  A,  B,  and  C  coefficients  found  by  least-squares  for  each  line  are 

shown  on  the  graph.  The  lines  were  fitted  to  the  strneture  values  only 

np  to  a  nean  separation  of  360  ha.  This  entoff  was  used  since  the 

strneture  beyond  this  point  decreased  and  would  require  the  fit  of  a 

■ore  coaplex  curve.  Other  curves  were  also  fit  to  these  points  as  will 

be  Mentioned,  and  the  cutoff  for  curve  fitting  with  diatance  was  a 

function  of  the  radiance  channel  and  synoptic  situation. 

The  linear  least-squeres  line  in  Figure  2.2  proved  to  have  a 

negative  structure  intercept  at  xero  separation  distance,  so  that  line 

was  not  considered  a  good  fit  since  the  aean-squared  difference 

(structure)  values  are  all  required  to  be  greater  than  xero.  However. 

the  sysmetric  quadratic  and  exponential  eurves  both  had  a  siailar 

2  —12 

intercept  of  0.12  (mV/(n  er  on  ))  it  0  km.  Ikie  translates  to  a  mean 

2  —1 

rms  noise  of  0.24  aW/(a  sr  on  )  sinoe  the  foraer  value  is  the  double 
aean-squared  noise*  one  ooaponent  coaing  froa  each  of  the  paried  values. 
This  value  is  considered  the  nost  probable  intercept  since  it  is 
confiraed  by  two  least-squares  lines.  A  potential  advantage  to  the 
syaaetric  quadratic  is  the  fact  that  its  slope  is  defined  as  xero  at 


13 


i? 

n 

Yj 


sero  separation  distance.  Die  structure  should  theoretically  decrease 
to  sneh  a  plateau  as  the  minimum  noise  level  is  reached.  However,  for 
soae  situations  the  plateau  aay  be  reached  quite  abruptly. 


Table  2.1  shows  the  results  for  each  of  the  two  regions  studied  as 


coapared  to  the  predicted  noise  levels  for  all  of  the  SSH-1  radiance 
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channels.  The  goal  here  was  to  see  how  the  obtained  radiance  noise 
levels  coapared  to  those  predicted.  The  two  regions  under  study  were 
kept  separate  to  see  if  eaoh  had  different  characterisitcs  or  if  one 
region  confined  the  values  obtained  by  the  other  region.  The  pressure 
at  the  weighting  funotion  aaziaua  for  eaoh  radianoe  channel  is  also 
shown.  Channels  which  peak  lower  in  the  ataosphere  are  ezpeeted  to  have 
larger  structure  functions  due  to  acre  ataospherie  and  surface 
variability.  This  ataospherie  variability  should  be  mostly  removed  by 
the  eztrapolation  process  to  zero  separation  distance.  However,  in  the 
case  of  lower  tropospheric  ohannels,  especially  the  window  channel,  this 
stmospheric  and  surface  variability  aay  not  be  coapletely  eliminated. 
Saall  cloud  aaounts  and  low  olouds  are  especially  hard  to  detect.  This 
will  be  addressed  further  in  a  following  section. 

The  noise  levels  obtained  by  structure  snalysis  coapared  favorably 
with  the  predicted  noise  levels  given  by  Yoder  (1975)  and  by  Klein,  et 
al  (1976).  The  nuabers  given  by  these  sourees  are  presumably  for  only 
one  satellite  and  the  results  here  are  composited  for  three  satellites. 
This  was  not  meant  to  be  an  ezhaustive  noise  analysis,  only  to  see  if 
the  noise  levels  obtained  here  were  grossly  different  than  the  predicted 
values,  which  did  not  seea  to  be  the  ease. 
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Hie  HjO  channel  noise  level*  for  the  aid-latitude  region  eoapared 

favorably  with  the  predieted  values.  Even  the  larger  predicted  noise 

for  the  F8  channel  (3S3  on  2)  did  show  up  in  the  structure  analysis. 

The  H^O  radiance  noise  levels  for  the  sub-tropic  region  were  all  about 
2  —1 

0.5  nW/(n  sr  cm  ).  This  behavior  is  unexplained  but  thought  to  be  a 
function  of  the  aore  moist  synoptic  situation  compared  to  the  aid- 
latitudes.  These  values*  however,  were  still  of  the  same  order  of 
magnitude  as  those  predicted. 

For  many  of  the  (X^  channels  the  noise  levels  for  the  mid-latitude 
and  sub-tropic  regions  showed  reasonable  comparison  to  predicted  values, 
however,  some  noise  levels  were  about  2  times  the  predicted  noise 
(reference  channels  E4  and  E5).  Many  of  these  noise  values  were 
obtained  by  aore  eoaplex  extrapolations  than  the  ones  for  the  example 
shown  in  Figure  2.2.  The  most  common  fit  was  for  the  form: 

Y  -  A  exp(Bx  +  Cx2) 

(2.4) 

This  is  not  a  symmetric  function  hut  it  did  many  times  provide  a  lower 
noise  level  than  the  other  functions  used  in  the  example.  A  ehoioe 
between  the  intercepts  for  all  the  fitted  curves  was  mad*  for  the  most 
probable  intercept  among  the  available  ohoices.  The  lowest  intercept 
was  oho  sen  if  any  one  value  was  not  confined  by  two  or  nor*  fitted 
lines. 

The  window  and  oxone  channels  also  showed  higher  noise  than 
predicted,  but  the  differences  are  not  unreasonable  since  the  window 
channel  especially  may  contain  some  unexplained  surfaee  variability  in 
adjacent  FOVs.  Only  gross  discrepancies  would  be  of  concern.  But, 
enlarged  noise  levels,  if  true,  may  signify  some  need  for  consideration. 


For  example,  •  minimum  information  retrieval  technique  depends  on 
statistical  knowledge  of  radianoe  noise  levels.  An  error  in  this  value 
does  effect  the  retrieval  solutions. 


2.4  Eigenvector/ value  Analysis 

Since  many  of  the  H^O  channels  are  rather  opaque  there  is  much 
redundancy  in  the  radiances  used  in  any  retrieval  program.  Therefore, 
an  analysis  of  the  same  limited  data  sets  as  used  in  the  structure 
analysis  was  performed,  which  showed  that  only  3  eigenvalues  explain 
over  99.5%  of  the  variance  in  the  IjO  channels.  The  other  5  eigenvalues 
add  only  an  insignificant  amount  of  information.  This  indicates  that  no 
more  than  3  independent  pieoes  of  information  exist  among  the  8  I^O 
channels. 

For  the  COj,  window,  and  oxone  channels  the  number  of  significant 
eigenvalues  needed  to  explain  over  99.5%  of  the  total  variance  is  again 
only  3.  This  means  that  no  more  than  3  independent  pieces  of 
information  exist  in  the  6  CO^,  window,  and  oxone  channels.  Much  of 
this  interdependence  among  ohannels  is  caused  by  the  tendency  for 
atmospheric  variables  to  be  correlated  in  the  vertical,  as  well  as  the 
broad  weighting  functions  for  each  channel. 


3.0  SCAN-ANGLE  CORRECTION  PROBLEM 

In  working  with  the  SSH-1  radiances,  the  problem  of  variations  due 
to  changing  satellite  viewing  angles  was  anticipated.  Some  radiance 
channels  would  experience  limb  darkening  (those  which  obtain  their  main 
contribution  from  the  troposphere)  and  others  would  experience  limb 
brightening  (those  which  obtain  their  main  contribution  from  the 
stratosphere).  This  introduces  an  extra  variation  in  addition  to 
changes  due  to  different  atmospheric  conditions  between  adjaoent 
satellite  measurements.  By  eliminating  this  variation  we  can  treat  all 
measurements  as  if  they  originated  at  a  common  (xero)  zenith  angle. 

The  software  documentation  provided  by  AFGWC  gave  a  3rd-order 
polynomial  correction  for  scan-angle  variations  of  the  form: 

Ratio  a  +  Aj  Zen  +  Aj  Zen^  +  A^  Zen^  (3.1) 

where  Ratio  is  the  factor  by  which  to  multiply  the  radiances  at  the 
specified  zenith  angle  to  correct  them  to  nadir.  Such  a  correction  is 
necessary  when  working  with  radianoes  from  radiometers  on  nearly  all 
satellites.  The  radiances  should  all  be  converted  to  a  common  (zero) 
zenith  angle,  or  alternately,  the  zenith  angle  variations  should  be 
taken  into  account  by  other  means,  such  as  in  the  retrieval  process. 

Typically  the  scan  angle  variations  follow  a  cosine  or  secant 
function  since  the  atmospheric  transmittances  are  a  function  of  the 
seoant  of  the  viewing  (or  zenith)  angle  with  respect  to  the  earth's 


normal  vector.  The  transmittance  calculations  follow  a  formula: 

x  »  t0  exp [see (Zen)]  (3.2) 

where  tq  is  the  transmittance  at  zero  zenith  angle  and  t  is  the 
transmittance  at  zenith  angle  Zen.  This  simple  relationship,  however, 
applies  only  to  a  single  atmospheric  level  and  a  more  oomplex 
relationship  evolves  as  the  radiative  transfer  equation  is  integrated 
through  the  different  density  layers  of  the  atmosphere. 

Past  work  of  the  authors  using  VTPR  data  revealed  that  the  zenith 
angle  correction  was  linearly  related  to  the  square  root  of  the  secant 
of  the  zenith  angle  (see  also  Duncan,  1977) .  The  secant  of  the  zenith 
angle  is  a  logical  choice  for  the  independent  variable  beoause  of  the 
above  explanation  relating  transmittance  variations  to  zenith  angle. 

For  the  SSH-1  radiance  channels  the  zenith  angle  corrections  as 
determined  by  polynomial  coefficients  provided  by  AFGWC  are  plotted  in 
Figure  3.1.  The  required  radianoe  changes  are  the  polynomial  ratios 
multiplied  by  the  mean  radianoe  in  eaeh  channel.  All  channels  except  F8 
(333  cm  *)  gave  a  near-zero  correction  at  zero  zenith  angle.  (There 
should  be  no  change  at  nadir;  such  a  change  would  bias  the 
measurements.)  However,  channel  F8  added  about  4  mV/(m  sr  cm  )  at  zero 
zenith  angle.  This  introduced  a  bias  of  about  +5%  to  the  mean  value  of 
about  73  mW/ (jz^*r  cm  ^)  for  this  channel. 

Even  with  this  bias  removed,  by  making  the  zero  order  polynomial 
coefficient  equal  to  1.,  the  F8  polynomial  had  a  shape  not  typical  of 
the  other  HjO  channels,  as  shown  in  Figure  3.2.  These  figures  were 
brought  to  the  attention  of  AFGWC  personnel.  Their  response  was  that 
the  polynomials  were  created  by  fitting  ourves  to  large  samples  of 
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satellite  measurements  averaged  at  each  soan  position.  The  polynomials 
were  based  on  data  from  satellite  F3  (not  to  be  confused  with  channel 
F3) .  The  average  F8  channel  response  with  zenith  angle  was  especially 
erratic  and  the  polynomial  which  was  fitted  to  the  zenith  angle  changes 
reflected  this  erratic  behavior.  No  other  explanation  for  the  F8  bias 
was  given. 

In  an  attempt  to  recreate  the  polynomials,  the  SSH-1  radiances  from 
at  least  two  time  periods  were  analzyed  to  determine  the  corrections 
necessary  to  take  the  zenith  angle  variations  out  of  the  several 
channels.  These  two  time  periods  were  6-7  November  1979  and  14-15 
January  1979.  The  analysis  required  considerable  computer  time  since 
each  period  consisted  of  many  orbits  of  data.  Some  problems  arose  in 
reading  the  November  1979  data,  but  the  January  1979  data  provided  a 
good  statistical  sample.  Over  55  thousand  aoan  spots  were  anlyxed  of 
which  over  27  thousand  were  considered  non-cloudy. 

Table  3.1  gives  the  mean  and  standard  deviation  for  each  of  the 
radiance  channels  as  compiled  from  the  14-15  January  1979  SSH-1  data 
supplied  to  us  by  AFGYC.  Effort  was  sude  to  eliminate  bad  radianoes 
[especially  those  caused  by  the  glare  obstructor  (glob)  problem 
(Valovcin,  1980)].  Also,  cloudy  radianoes  were  eliminated  by  using  a 
minimum  threshold  of  85  mW/(m  sr  cm  )  for  the  12  pm  window  channel 
(835  cm  *) .  Radiance  values  lower  than  this  threshold  were  considered 
to  be  cloud  contaminated. 

After  the  cloudy  values  were  eliminated  by  the  threshold  method, 
polynomials  similar  to  Equation  3.1  were  fit  to  this  data  and  are 
plotted  in  Figure  3.3.  The  required  radiance  changes  are  similar  to 
those  provided  by  AFGYC.  However,  most  of  the  F  channels  show  similar 


Table  3.1 


SSH-1  radiaaoe  etatiaties  for  14-15  Jaa  1979  aad  eorraetioaa  at  S7° 
sealtb  angle  for  polyaoaials  fit  to  tke  Jaa  79  data  sat  aad  for 
polyaoaials  prorlded  by  AF6VC 

1*1  19  PolraoMlal  Correction  at  37° 


Channel 

Vavenaaber 

Msaa 

Std.  Dot. 

Jaa  79 

AFOVC 

(oa-1) 

<a»/(a2 

sr  ca”1)) 

F8 

353 

73. 

4.0 

3.4 

5.2 

FI 

355 

77. 

3.9 

3.3 

3.6 

F7 

374 

84. 

4.4 

3.6 

4.1 

F2 

397 

88. 

4.6 

3.7 

4.7 

F< 

408 

99. 

5.3 

3.7 

5.3 

F3 

420 

95. 

5.2 

3.7 

5.1 

F4 

441 

101.  ' 

5.7 

4.1 

5.7 

F5 

535 

116. 

6.6 

4.2 

6.3 

El 

6(8 

57. 

3.0 

-2.3 

-2.1 

E2 

677 

44. 

2.6 

-2.2 

-1.5 

E3 

695 

42. 

2.6 

1.4 

0.3 

E4 

708 

58. 

2.9 

7.7 

6.5 

E3 

725 

78. 

4.6 

9.6 

9.7 

E 6 

747 

91. 

5.9 

8.7 

10.4 

VI 

835 

(two  looks) 

108. 

10.7 

3.5 

— 

V2 

835 

(oae  look) 

108. 

10.7 

3.4 

6.6 

change*  and  ara  clustered  together  with  the  window  channels  VI  and  W2 
(two  and  one  looks  respectively).  The  F  channel  xenith  angle  ohanges  in 
Figures  3.1  and  3.2  are  aare  spread  out.  bnt  thia  difference  in  spread 
is  aost  likely  due  to  the  different  data  samples  which  went  into  eaoh 
analysis. 

The  last  two  eoluams  of  Table  3.1  give  the  polynomial  corrections 
specified  for  the  maximum  xenith  angle  of  57°.  The  'Jan  79'  column 
gives  the  correction  for  polynomials  fit  to  the  SSH-1  data  for  which  the 
statistics  were  given  in  the  third  and  forth  oolnans  of  Table  3.1.  The 
'AFGVC'  column  gives  the  same  correction  (at  57°)  for  the  polynomials 
provided  by  AFGWC  (on  an  unknown  data  set).  The  'Jan  79'  data  inoludes 
passes  from  two  satellites  (VX3536  and  VX4537),  whereas  the  'AFGVC' 
polynomials  were  fit  to  only  one  satellite  (F3) .  Regardless  of  these 
differences,  the  corrections  specified  by  the  polynomials  for  the  two 
data  sets  are  within  about  2  aV/(a?sr  cm-1)  for  most  of  the  F  (^0)  and 
E  (C02)  channels.  The  largest  disorepanoies  occurred  for  the  most 
transparent  channels  (F5,  B6,  and  V2)  and  for  ohannel  F8  where  the  AFGVC 
polynomial  is  assumed  to  be  in  error. 

A  solution  to  this  problem  is  to  assume  that  the  polynomial 
correction  for  F8  is  similar  to  that  for  FI,  the  nearest  H2<)  channel  in 
terms  of  wavenumber  and  sensing  depth  in  the  atmosphere.  The  required 
radiances  changes  for  the  Jan  1979  data  are  similar  for  channels  F8  and 
FI,  but  rather  than  using  the  polynomials  fit  to  this  data  set  it  was 
assumed  to  be  safe  to  apply  the  AFGVC  polynomial  for  FI  to  channel  F8. 
This  produces  a  radiance  correction  of  3.4  (instead  of  the  previous 
value  of  5.2)  mV/(a2sr  cm  *)  at  the  maximum  xenith  angle  of  57°. 

It  was  also  decided  that  the  polynomials  used  to  oorrect  the  SSB-1 


radiances  should  produce  no  bias  at  aero  senith  angle.  Some  snail 
biases  were  produoed  by  the  APS VC  polynonials.  This  bias  was  elininated 
by  always  asking  the  aero  order  ooeffieient  equal  to  1.  The  new 
polynonials  (with  no  bias)  and  the  eorreoted  channel  F8  polynonial 
(equal  to  FI)  are  shown  in  Figure  3.4.  The  required  corrections  for 
channels  F8  and  FI  show  sons  slight  differences  ewen  though  the 
polynonials  (ratios)  are  identical.  The  absolute  correction  is 
dependent  on  the  nean  walues  used  to  plot  Figure  3.4.  The  aean  value 
for  channel  FI  is  larger  than  the  aean  for  channel  iT8  (as  shown  in  Table 
3.1).  therefore  the  required  correction  is  sinilar  but  not  identical  for 
each  channel. 

These  zenith  angle  correction  polynonials  are  now  better  suited  for 

application  without  biasing  or  adversely  affecting  SSB-1  statistics. 

However,  it  is  suggested  that  any  senith  angle  correction  be  aore  data 

sanple  speoific  or  be  general  enough  to  apply  latitudinally  and/or 

seasonally.  The  differences  between  the  required  corrections  at  57°  for 

the  two  data  sets  for  soae  of  the  acre  transparent  ohannels  is  auch 

larger  than  the  radiance  noise  level  (assuring  a  noise  level  of  about 
2  —I 

0.23  aV/(n  sr  ca  )).  However,  the  polynoaials  provided  by  AFOVC.  with 
the  improvements  suggested  here,  do  produce  a  required  correction  for 
channel  F8  that  does  agree  with  the  January  1979  data  saaple  studied. 

The  suggestions  aade  here  will  be  applicable  to  SSH-2  data  which  is 
to  be  available  in  1982.  These  zenith  angle  noaalizing  procedures  will 
continue  to  be  required  if  a  statistical  retrieval  scheme  is  to  be 
applied  uniformly  to  radiance  measurements  from  various  viewing  angles. 
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Figure  3.4.  Required  SSH-1  radianoe  change  as  a  function  of  xenith 

angle  for  AFGVC  polynomials  with  channel  F8  polynomial  equivalent 
to  that  of  channel  FI  and  for  no  bias  at  xero  xenith  angle. 
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4.0  CLOUD  ELIMINATION  TECHNIQUES 


la  both  the  structure  analysis  and  the  scan-angle  correction,  it  is 
important  to  eliminate  the  effects  of  clouds  on  the  SSH-1  radianoes. 
Clouds  normally  reduce  the  radiance  by  an  amount  which  is  proportional 
to  the  cloud  amount,  ao  varying  amounts  of  cloud  in  adjacent  FOVs  cause 
the  structure  to  be  abnormally  high.  Tith  clouds  eliminated  the 
etrueture  should  reflect  only  mean  horizontal  gradient  in  the 
atmospheric  and  aurface  parameters. 

One  way  to  eliminate  clouds  is  to  establish  a  minimum  threshold 
radiance  below  which  the  FOV  is  assumed  to  contain  at  least  a  fraction 
of  cloud.  Usually  this  threshold  is  picked  to  correspond  to  a  certain 
surface  temperature.  The  window  channel  at  11  pm  waa  used  here  since  it 
is  the  most  transparent  to  the  atmosphere  and  most  affected  by  clouds. 

In  the  case  of  the  mid-latitude  area  under  study,  thresholds  of  75,  80, 
and  85  rif/(m*sr  cm  )  were  chosen.  These  correspond  to  surface 
temperatures  of  approximately  -8°,  -4°,  and  0°C.  These  thresholds  were 
tested  experimentally  and  the  exact  value  used  should  depend  on  the 
expected  background  (surface)  temperature.  For  the  mid-latitude  region 
with  measurements  taken  in  November  the  lowest  threshold,  corresponding 
to  -8°C,  seems  to  be  optimal,  whereas  much  higher  threshold  temperatures 
were  used  for  the  subtropic  region  (typically  0  or  10°C) .  Surfaces 
colder  than  the  threshold  probably  contain  at  least  a  small  fraction  of 
cloud,  causing  the  observed  radianoe  to  correspond  to  a  lower  equivalent 
temperature  than  the  actual  surfaee  temperature. 


A  second  cloud  detection  approach  vas  also  nsed.  This  method 
relied  on  the  two  measurements  made  using  the  window  channel  filter. 
(Each  SSH-1  radiance  channel  is  sensed  twice  for  each  FOV) .  This  method 
is  a  variation  of  the  technique  to  determine  surface  temperature  and 
cloud  amount  using  two  window  ohannels  in  different  spectral  regions 
(Smith  and  Rao,  1973),  which  will  be  discussed  in  the  next  subsection. 

Using  the  two  looks  through  the  window  filter  it  was  thought  that 
any  slight  shift  in  the  FOV  or  change  in  cloud  characteristics  (if 
possible)  would  be  reflected  in  a  ohanging  window  radiance.  The  two 
looks  would  give  slightly  different  radiances  since  the  observed  cloud 
field  for  one  look  would  be  slightly  different  than  for  the  other  look. 
This  method  was  used  to  eliminate  FOVs  with  normally  high  radiances, 
above  the  minimum  threshold,  but  which  possibly  contained  small  cloud 
amounts.  Ihe  only  conf Iruatioa,  however,  was  that  this  method  seemed  to 
eliminate  those  spots  whioh  had  lower  (and  possibly  cloudy)  radiances 
when  compared  to  adjacent  scan  spots.  The  structure  values  became  more 
well  behaved  once  this  detection  method  was  in  effect,  in  that  much 
small  scale  variability  seemed  to  be  eliminated,  and  the  curves  fit  to 
the  1-dimensional  structure  values  produced  more  realistic  noise  levels. 

In  the  case  of  the  two  window  looks,  the  radiance  for  one  look  is 

given  in  data  provided  by  AFGVC,  and  the  mean  radiance  for  the  two  looks 

is  also  given,  so  the  difference  between  one  look  and  the  mean  would  be 

half  that  between  the  two  measurements.  This  is  just  mentioned  because 

it  is  thought  it  would  be  better  to  have  both  measured  values  rather 

than  one  value  and  the  mean  sinoe  the  difference  is  normally  not  very 
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large,  usually  on  the  order  of  0.5  mV/(m  sr  cm  )  or  less.  Radiance 

2  -1 

difference  thresholds  of  0.1  or  0.2  mW/(m  sr  cm  ),  therefore,  seemed  to 


eliminate  most  small  scale  noise 


A  number  of  different  combinations  of  the  two  thresholds  (minimum 
value  and  difference)  were  used  in  the  structure  analysis  provided  in 
the  previous  sections.  For  ^btple,  the  structure  values  in  Figure  2.2 
were  calculated  using  window  channel  thresholds  of  100  *W/(m*sr  cm-*) 
minimum  and  0.1  mW/(m  sr  cm  )  difference.  The  noise  level  results 
(extrapolations)  were  dependent  on  the  thresholds  chosen.  A  number  of 
noise  levels  were  then  generated  and  the  most  probable  value  was  chosen 
when  that  value  was  obtained  using  different  threshold  combinations. 
Usually  the  tighter  thresholds  (higher  minimum  window  radiances  and 
lower  window  difference)  produced  more  reasonable  noise  results.  The 
lower  noise  level  was  chosen  in  oases  where  an  alternate  but  higher 
noise  level  seemed  unreasonable. 

4.1  Slnalo  FOV  Cloud  Determination  Psina  Two  Window  fflUMMUll 

By  using  coincident  measurements  from  two  window  channels  with 
different  spectral  responses  we  can  effectively  split  a  single  Field- 
of-View  (FOV)  into  a  dear  sub-area  and  an  adjacent  cloudy  sub-area  each 
with  a  different  temperature.  This  capability  is  due  to  the  different 
spectral  responses  of  each  channel.  Various  assumptions  must  be  made 
about  the  FOV  and  the  window  channels.  The  first  assumption  is  that 
only  two  separate  sub-areas  comprise  the  total  FOV  (i.e.  the  surface  and 
cloud  sub-areas  are  uniform).  A  second  assumption  is  that  the  differing 
measurements  between  window  channels  are  more  a  result  of  special 
variations  in  the  viewed  temperature  field  than  of  any  atmospheric 
absorption  differences  between  the  two  window  channel  values.  The 
equation  governing  either  of  these  window  ohannels  is  then: 


where  B  is  the  Plaaek  function,  k^  is  the  channel  wavenumber,  and  the 
FOV  is  split  into  a  clear  and  cloudy  fraotion  with  p  equal  to  the  cloud 
amount.  As  their  name  isqplies,  these  window  channels  encounter  little 
atmospheric  absorption.  Therefore,  the  two  terms  in  the  equation 
originate  from  the  surface  and  cloud  fractions,  respectively,  with  no 
atmospheric  contribution. 

If  two  window  channels  such  as  3.7  and  11  |im  are  used  to  view  such 
a  split  FOV,  then  the  individual  speotral  responses  of  the  channels  will 
cause  eaoh  to  give  a  different  effective  temperature,  ?eff  This  is  due 
to  the  non-linear  resonse  of  the  Plank  function  with  respect  to 
temperature  variations.  On  the  other  hand,  a  uniform  FOV  would  give  the 
same  effective  temperature  in  both  channels.  For  example,  a  uniform  0°C 
(273  K)  FOV  will  give  radiances  of  0.18  and  76.3  mW/(m2sr  cm-1)  at  3.7 
and  11.1  pm,  respectively.  However,  a  FOV  whioh  is  partially  cloud 
contaminated  would  give  lower  radiances  and  lower  effective  temperatures 
in  both  channels.  The  channel  affected  most  would  be  the  longwave  (11.1 
pm)  channel.  The  warm  (surfaoe)  fraotion  of  the  FOV  would  effectively 
contribute  more  to  the  3 .7  pm  radiance  than  the  11 .1  pm  radianee 
although  both  channels  see  the  same  FOV  and  therefore  the  same  cloud 
amount. 

The  use  of  the  two  window  ohannels  for  cloud  determination  is  not 
without  precedence.  Smith  and  Rao  (1973)  dealt  with  the  3.7  and  11  pm 
window  channels  as  applied  to  two  adjacent  FOVs.  They  solve  four 
equations  with  four  unknowns:  the  surfaoe  temperature,  the  cloud 
temperature,  and  the  oloud  amount  in  each  FOV.  To  solve  this  set  of 


equations  they  a sane  that  the  average  surface  and  cloud  teaperatures  are 
the  saae  for  the  tvo  FOVs. 

Dozier  (1980),  and  Matson  and  Dozier  (1981),  however,  dealt  with 
the  saae  window  channels  using  only  one  FOV  by  assuming  that  we  know  or 
can  assume  the  surface  teaperature,  T  This  requires  only  two 
equations  to  be  solved  for  the  cloud  temperature  T  ^  and  cloud  amount 
p.  The  solution  is  dependent  on  the  assumed  surface  temperature,  but 
the  solution  is  unique.  The  noise  on  the  radiances,  however,  may  place 
the  solution  outside  of  physically  reasonable  limits. 

By  eliminating  the  cloud  amount  p  from  the  two  Equations  4.1  we 

get: 

A1  '  *  v  »«cia.kj)  ♦  A3  -  0  (4.2) 


where  Aj  and  A3  are  known  if  the  surface  teaperature  is  assumed. 


*1  “  Bsfc,2  “  3eff,2 
A2  “  Beff,l  _  Bsfc,l 
A3  "  Beff,2  *  Bsfc,l  "  Beff»l  Bsfc,2 

and 


Bsfc, i 
Beff,i 


*  B<T.fc'ki> 

-  B<Teff.i'V 


The  effecitve  temperature  T^^  in  each  channel  is  a  known  measurement, 
as  well  as  the  channel  wavenumber  k^. 

Equation  4.2  cannot  be  solved  explicitly  since  it  is  nonlinear  in 


temperature.  An  iterative  solution  is  siaple  if  a  starting  cloud 


temperature  it  assumed.  A  simple  starting  assumption  would  be 


old 


<  T,.  <  T. 


lll 


3.7 


<  T 


sfc 


(4.3) 


After  the  cloud  temperature  is  found*  the  cloud  fraction  is 
determined  by  solving  for  p  from  the  original  two  equations,  one  for 
each  channel 


where 

A4  =  Bcld,2  *  Bsfc,l  ~  Bcld,l  *  ®sfc,2 


(4.5) 


and 


Bcld, i 


B(Tdd'ki> 


(4.6) 


If  there  were  equal  signs  in  Equation  4.3  the  clond  amount  p  either 
equals  0  or  1.  These  limiting  oases  yield  a  result  involving  no 
distinction  between  the  two  window  channels  (i.e.  a  FOV  comprised  of 
only  a  single  effective  temperature,  either  totally  clear  or  totally 
cloudy) . 

The  assumption  of  little  or  no  atmospheric  contribution  for  the 
window  channels  can  be  overcome  by  correcting  the  window  channels  for 
any  atmospheric  absorption  which  might  ooour.  This  can  be  done  by  using 

i 

assumed  atmospheric  temperature  and  moisture  profiles  to  estimate  the 
atmospheric  transmittance  in  each  channel. 

!  The  radiative  transfer  equation  (RTE) 

I  *r. 

!  UT^)  -  B(Tef£,kA)  rl(T(p),Q(p))  +  X  BtKp),^)  dp 

! 


(4.7) 


can  be  solved  for  the  effective  temperature  T  _  .  given  the  equivalent 
brightness  temperature  T^.  The  transmittance  xi  for  each  channel  is 
calculated  from  the  assumed  temperature  and  moisture  profiles  T(p)t 
Q(p).  These  atmospheric  corrections  are  small  inmost  cases,  typically 
1  to  2  K  for  a  relatively  dry  atmosphere. 

Another  necessary  correction  would  be  for  reflected  radiation  in 
the  shortwave  (3.7  pm)  channel.  This  correction  can  be  on  the  order  of 
10  K  for  50%  cloud  cover.  (See  Hayden,  et  al  (1981)  for  details).  This 
reflected  correction  should  be  applied  before  any  atmospheric  correction 
is  made. 

This  single  FOV  technique  has  been  tested  on  some  TIROS  Operational 
Vertical  Sounder  (TOYS)  radiances  from  the  TIROS-N  and  NOAA-6 
satellites.  Figure  4.1  shows  the  window  brightness  temperature 
difference  as  a  function  of  the  satellite-derived  bi-directional 
reflectance.  The  high  correlation  of  0.89  for  this  data  signals  a 
strong  relationship  and  the  ability  of  the  two  window  channels  to  detect 
clouds  based  on  their  difference.  The  same  techniques  can  be  applied  to 
the  SSB-2  instrument  which  will  contain  both  3.7  and  11  |im  window 
channels.  Cloud  temperature  (or  cloud  height)  and  cloud  amount  can 
therefore  be  determined  for  a  single  FOV.  Uncertainty  in  the  radiances 
and  the  corrections  to  the  radiances  will  cause  some  non-physical 
solutions  (i.e.  p  <  0  or  p  <  1) .  These,  however,  can  be  handled  as 
limiting  cases  p  ■  0  or  p  -  1  within  the  bounds  of  the  uncertainty  in 


the  measurements  and  assumptions 
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Figure  4.1.  Hie  window  brightness  temperature  difference  (WBTD)  as  a 
function  of  the  satellite-derived  bi-directional  reflectance  for 
NOAA-6  HISS-2  window  channels  at  3.7  end  11  pa.  A  correlation 
coefficient  of  0.89  is  given. 


5.0  STATISTICAL  RETRIEVAL  ASSESSMENT 

The  sain  thrust  here  was  to  test  s  ststistiesl  eigenvector 
retrieval  technique  such  as  used  at  AFGVC  for  SSH-1  retrievals.  This 
algorithm  was  compared  to  simple  regression  to  deteraine  the  maximum 
retrieval  capabilities  on  the  dependent  data  set.  An  independent  (test) 
data  set  would  give  larger  'errors'.  The  'errors*  here  are  differences 
between  RAOB  temperature  and  aoisture  values  and  those  retrieved  at 
nearby  satellite  soundings.  It  should  also  be  noted  that  the  'errors' 
mentioned  in  this  section  include  errors  in  the  RAOBs  and  differences  in 
the  time  and  spatial  coverage  of  satellite  soundings  versus  RAOBs. 
Differences  in  tine  of  up  to  +  3  hours  were  allowed,  as  were  differences 
of  up  to  100  nautical  ailes  (185  ha)  in  location  of  the  RAOB  versus  the 
soan  spot  center.  The  statistical  data  base  aoquired  from  AFGVC 
contained  542  such  pairs  which  eomposed  the  data  set  used  in  this  studty. 
All  regions  of  the  globe  are  covered,  but  the  soundings  are  weighted 
heavily  towards  the  northern  heaisphere. 

In  an  effort  to  deteraine  the  Mxiaua  accuracy  at  which 
teaperatures  and  nixing  ratios  can  be  retrieved  from  the  SSH-1  data,  two 
computer  programs  were  written  to  do  statistical  retrievals.  The  first 
does  retrievals  by  multiple  linear  regression.  It  works  as  follows. 

Let  jt  be  a  vector  of  15  temperature  deviations  from  the  mean  at  the  15 
standard  levels  from  10  to  100  Ub  and  tg  be  a  vector  of  15  brightness 
temperature  deviations  from  the  mean  for  the  6  CO^,  8  H^O,  and  one 
window  channel  on  the  SSH-1.  Then  we  want  to  find  a  matrix  C  which 


relates  the  two  in  the  Benner: 


~  “  %  (5.1) 

The  aetrix  C  is  found  by  the  least- squares  method*  that  ia*  by  forming 
the  aatriees  of  T  and  whose  eolnans  are  the  £  and  vectors 

respectively  from  the  paired  radianee-RAOB  data.  Ihe  best  fit  C  is  then 

C  -  TlJCTglJr1  (5.2) 

where  the  superscript  T  indicates  transpose*  and  the  superscript  -1 
indicates  the  inverse.  Regression  of  airing  ratio  was  done  similarly* 
but  with  values  of  airing  ratio  in  £  vectors  of  length  (. 

For  the  first  experiaeat*  matrix  C  was  calculated  for  data  between 
30°N  and  60°N  (aid-latitude)  and  far  dates  between  26  June  1979  and  2 
February  1980,  which  included  118  of  the  542  satellite-RAOB  pairs.  (No 
desoriaination  was  aade  for  the  data  froa  different  SSB-1  instruments  on 
different  satellites).  Errors  in  the  regression  equation  using  this 
dependent  data  were  then  calculated.  This  error  is  an  estiittte  of  the 
best  that  the  SSH-1  instrument  is  able  to  do  at  retrieving  ataospherie 
paraaeters.  Figure  5.1  shows  the  teaperature  errors*  and  Figure  5.2 
shows  the  aixing  ratio  errors.  In  the  process  of  calculating  the  mixing 
ratios  some  extremely  large  values  were  found.  Figure  5.2  was 
calculated  after  aixing  ratios  greater  than  the  aean  plus  three  standard 
deviations  were  removed  froa  the  saaple.  These  large  aixing  ratios  nay 
or  aay  not  be  in  error,  but  in  estiaating  the  best  that  the  sounder 
could  do  it  seeaed  appropriate  to  throw  out  the  large  data. 

As  can  be  seen  in  Figure  5.1  an  S8B-1  instrument  oan  do  a  fairly 
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Figure  5.2.  Root-aean-square  airing  ratio  error  ee  a  function  of  pres- 
aure  for  the  SSH-1  regression  retrieval  on  the  saae  (dependent) 
data  set  used  to  develop  the  regression  coefficients. 


good  job  of  estinating  tenperatnres.  Correlation  coefficient*  between 
retrieved  tenperatnres  and  RAOB  teaperatnrea  are  0.9  or  above  ezoept 
near  the  tropopanae  (200  -  300  nb) .  Abaolnte  error*  are  lea*  than  3  K 
m  in  the  niddle  troposphere  (700  -  300  nb)  and  in  the  stratosphere  (10 
-  150  nb) .  Correlation  ooeff ieients  between  retrieved  and  observed 
nixing  ratios  are  not  as  high  as  for  tenperatnres.  and  the  error  at  low 
levels  is  1  -  2  g/kg.  If  this  technique  were  tested  on  data  other  than 
that  nsed  to  develop  the  regreasion  coeffioienta  (independent  data),  the 
results  would  be  worse  overall,  and  we  auspeet  they  would  be 
signif ieantly  worse  for  the  nixing  ratios. 

The  second  retrieval  progran  tested  does  retrievals  by  the  Snith- 
Voolf  eigenvector  technique  (Snith  and  Voolf,  1976) .  The  purpose  of 
this  was  to  find  any  errors  in  the  eigenvector  technique  whioh  is  in 
operational  uae  at  AFGfC.  If  the  nunber  of  eigenvector*  uaed  in  the 
solution  is  equal  to  the  nunber  of  brightness  tenperatnres  (channels) 
then  the  eigenvector  technique  is  identical  to  the  regression  technique. 
If  fewer  eigenvectors  are  used  to  deternine  the  C  natrix,  the  error  is 
increased  for  dependent  data.  (The  eigenvector  technique  oan  actually 
outperfom  the  regression  technique  on  independent  data).  Shown  in 
Figures  5.3  and  5.4  are  the  nass-weighted  average  rns  errors  in 
retrieving  teaiperatures  and  nixing  ratios,  respectively,  versus  the 
nunber  of  eigenvectors  used  to  deternine  the  C  natrix.  Shown  for 
eonpariaon  is  the  nass-weighted  average  rns  error  for  the  regression 
technique. 

The  naxinun  possible  nunber  of  eigenvectors  whioh  can  be  used  to 
construct  the  C  natrix  is  the  lessor  of  the  nunber  of  levels  to  be 
retrieved  and  the  nunber  of  channels  used  in  the  retrieval  —  15  for 
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Figure  5.1.  Eoot-aeaarsqnare  tenperatnre  arror  aa  a  function  of  prea- 
anra  for  the  SSH-1  regression  retrieval  on  the  sane  (dependant) 
data  sat  need  to  develop  tha  ragraaslon  coefficients. 
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Figure  5.4.  Mass-veighted  average  ns  airing  rstio  error  ee  e  function 
of  the  nnaber  of  eigenvectors  need  in  the  SSH-1  eigenvector 
retrieval. 


temperatures,  6  for  mixing  ratios.  Figures  5.3  and  5.4  show  that  using 
about  half  of  the  maximum  number  of  eigenvectors  (7  for  temperatures,  3 
for  mixing  ratios)  reduees  the  'error*  by  90%  of  the  potential 
reduotion. 

These  moisture  results  oompare  favorably  with  previous  tests  of  the 
accuracy  of  retrievals  using  the  SSH-1  radianoes  (Klein  et  al,  1976; 
Vachtmann  et  al,  1974;  Mount  et  al,  1977)  which  show  ms  temperature 
'errors'  of  2  to  4  K  as  typical.  However,  the  larger  'errors'  in 
retrieval  of  moisture  parameters  at  AFGVC  are  most  likely  not  caused  by 
the  retrieval  programs  used,  but  are  probably  a  function  of  the 
statistical  radiance-RAOB  sample  or  of  the  independent  data  sets  under 
study.  Some  large  moisture  variances  were  found  in  the  RAOB  data  and 
these  could  adversely  affect  the  retrieval  results. 

If  statistical  retrieval  schemes  continue  to  be  used,  then  the 
radiance-RAOB  pairing  should  be  better  scrutinised  to  eliminate 
suspected  noise  in  both  the  radiances  and  the  conventional  soundings. 
This  is  critical  to  the  statistical  retrieval  seheme  in  both  the 
production  and  use  of  the  C  matrix.  Casnsunication  with  personnel  at 
AFGVC  confirmed  that  suoh  problems  have  existed  and  should  be  corrected 
if  reasonable  moisture  results  are  to  be  obtained.  It  is  also  suggested 
that  the  temperature  and  moisture  retrievals  be  conducted  separately  in 
order  to  avoid  contaminating  the  retrieved  temperatures  with  noise  from 
the  moisture  channel  radianoes  and  vice  versa. 
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Regression  retrievals  were  also  performed  on  independent  data  sets 
(data  sets  other  than  those  used  to  develop  the  regression  equations). 
The  errors  for  retrievals  using  independent  data  should  inorease  above 


those  for  retrievals  using  dependent  data  sets.  In  this  application 
alternate  retrieval  schemes  can  nake  improvements  over  simple  regression 
retrievals.  An  eigenvalue/vector  retrieval  scheme,  for  example,  can 
actually  ont-perform  the  regression  retrievals  on  independent  data. 

This  is  done  by  eliminating  some  highr-order  eigenvectors  which  carry 
only  noise  information. 

The  independent- set  regression  experiments  were  oondncted  using  all 
542  souning  pairs  provided  by  AFGWC.  To  determine  the  minimum  error 
levels  obtainable  with  this  larger  set  of  data,  the  satellite-RAOB  pairs 
were  first  used  both  to  develop  the  regression  equations  and  then  as  a 
dependent  set  to  test  the  regression  retrieval  scheme.  Alternately,  to 
provide  both  a  dependent  and  an  independent  data  set,  the  542 
satellite-KAOB  pairs  were  split  various  ways  into  two  equal  size  sets. 
Four  methods  were  used  to  make  this  division.  Each  case  allowed  half  of 
the  of  the  set  (271  satellite-RAOB  pairs)  to  be  the  dependent  set  for 
development  of  the  regression  coefficients  and  the  other  half  to  be  the 
independent  set  for  testing  of  the  regression  equations.  The  divisions 
were : 

1)  odd  pairs  as  dependent  set,  even  pairs  as  independent  set; 

2)  even  pairs  as  dependent  set,  odd  pairs  as  independent  set; 

3)  first  half  as  dependent  set,  second  half  as  independent  set; 

4)  second  half  as  dependent  set,  first  half  as  independent  set. 

The  results  using  each  of  the  4  divisions  of  the  satellite-RAOB 
data  were  then  averaged  to  give  representative  'errors'  for  a  regression 
retrieval  scheme  on  data  other  than  that  used  to  derive  the  regression 
coefficients.  The  results  are  given  in  Tables  5.1  and  5.2  which  show 
the  temperature  and  mixing  ratio  errors  respectively.  Some  of  the 
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tabulated  results  are  also  plotted  in  Figures  5.5  and  5.6.  The  standard 
deviations  (square  root  of  total  variance)  of  the  aeteorologieal 
parameters  at  various  levels  are  also  given.  These  standard  deviations 
show  the  variability  of  the  data  and  oan  be  used  as  a  baaeline  for 
understanding  the  noise  levels.  The  explained  variance  values: 

E*"1’  V*r-  -  1  -  Uttl  vtriapo*  <5-3> 

give  a  measure  of  the  amount  of  variance  of  the  whole  data  set  that  can 
be  explained  by  the  regression  retrieval  scheme.  This  value  will  range 
from  0  to  1,  with  zero  indicating  that  the  retrievals  explain  none  of 
the  data  variability  (i.e.  the  error  variance  is  as  large  as  the 
variance  of  the  total  data  set).  An  explained  varianoe  value 
approaching  1  would  indioate  an  error  varianoe  much  smaller  than  the 
variance  of  the  total  data  set,  or  good  satellite  retrievals  when 
compared  to  RAOBs. 

The  results  in  the  first  four  oolumns  of  Table  5.1  show  that  the 
temperature  errors  reach  a  minimum  of  about  2  K  rms  in  the  lower 
stratosphere.  The  sample  standard  deviation,  however,  is  also  low  here, 
but  the  explained  variance  is  still  very  high.  Temperature  errors 
increase  to  4  I  rms  at  the  surface  and  also  increase  at  higher  levels. 
The  minimum  explained  varianoe  of  70%  occurs  at  the  25  kPa  (250  mb) 
level,  near  the  tropopause. 

The  error  statistics  for  the  various  half  dependent-half 
independent  data  set  divisions  are  also  shown  as  well  as  the  half 
dependent-half  independent  mean  errors.  These  errors  increase  slightly 
over  the  totally  dependent  sample  errors  and  range  from  2  K  in  the 
stratosphere  to  4.7  K  at  30  kPa,  with  errors  above  4  K  down  to  the 
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Figure  5.5.  Root-*® in- square  temperature  error  as  a  function  of  pres¬ 
sure  for  SSH-1  regression  retrievals  on  both  dependent  and  indepen¬ 
dent  (test)  data  sets. 
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Figure  5.6.  Root-seen- square  airing  ratio  error  as  a  function  of  pres¬ 
sure  for  SSH-1  regression  retrievals  on  both  dependent  and  indepen 
dent  (test)  data  seta. 


surface.  These  larger  errors,  however,  still  explain  at  least  55%  of 
total  varianoe  at  30  kPa  and  larger  anounts  at  other  levels. 

These  temperature  error  reaults  compare  favorably  with  previous 
tests  of  the  aceuraey  of  retrievals  using  the  SSH  radiances  (Mount,  et 
al.,  1977;  Klein,  et  si.,  1975)  which  show  ras  errors  of  2  to  4  K  as 
typical.  In  order  to  prove  aore  effective  than  simple  regression,  any 
retrieval  scheae  developed  will  have  to  provide  reduced  'errors'  in  a 
similar  satellite-KAOB  coaparisons. 

The  mixing  ratio  errors  for  the  6  standard  levels  are  shown  in 
Table  5.2,  in  a  foraat  aiailar  to  Table  5.1.  Because  aoisture  resides 
most  heavily  in  the  lowest  layers  of  the  atmosphere  the  'errors*  can  be 
expected  to  decrease  with  height.  The  explained  variance,  therefore, 
serves  especially  well  to  coapare  the  'error'  to  the  saaple  standard 
deviation  at  each  level,  putting  the  'error'  in  perspective.  Coluans  3 
and  4  give  the  all  dependent  saaple  results,  where  it  can  be  seen  that 
nixing  ratio  errors  increase  to  1.34  g/kg  at  the  surface.  Explained 
variances  at  all  levels  are  above  45%  with  the  largest  value  of  77%  at 
the  surface.  For  the  half  dependent-half  independent  case  the  aean 
errors  increase  to  1.48  g/kg  at  the  surface.  Explained  variances  here 
range  from  31  to  72%,  again  with  the  largest  percentage  explained 
variance  at  the  surfaoe.  This  is  nost  likely  due  to  the  larger  aoisture 
signal-to-noise  ratio  at  the  surfaoe.  While  the  nixing  ratio  error 
decreases  slightly  with  height,  it  does  not  deorease  with  height  nearly 
as  rapidly  as  the  saaple  standard  deviation. 

The  aoisture  results  in  general  show  a  reasonable  but  United 
aoisture  retrieval  capability.  Explained  variances  below  50%  show  very 
little  usable  information.  The  best  aoisture  information  therefore 


comes  from  n ear  the  surface.  Similar  explained  variance a  were 
determined  by  Savage  (1980)  Being  actual  the  AFGVC  atatiatieal 
eigenvector  retrieval  teehnique.  Other  moiatnre  retnlte  by  Crane  (1976) 
provide  explained  variance a  in  total  preeipi table  wager  of  elightly 
better  than  50%  ueing  the  H^O  channel a  only.  The  ability  to  derive 
moiatnre  information  vaa  alao  ah  own  to  be  enhanced  by  maing  aotual 
temperature  information  aa  independent  regreaaion  prediotora.  The 
explained  varianoea  then  roae  to  the  over  80%  range.  The  nae  of  COj 
channel  radianoea  direotly  inatead  of  aetttal  temper a tore a  would  degrade 
theae  reaulta  aomewhat. 

The  reaulta  uaing  thia  ainple  regreaaion  retrieval  scheme, 
therefore,  do  confirm  the  reaulta  of  aimilar  atudiea  by  othera.  The 
reaulta  of  courae  depend  on  the  data  aample  under  study,  but  with  large 
data  aamplea  the  reaulta  ahould  be  largely  invariant.  Any  improvement a 
to  be  made  could  poaaibly  ariae  out  of  a  more  phyaioal  retrieval  eehene 
whioh  would  better  interpret  the  atrengthe  of  eaeh  radianoe  channel.  An 
attempt  at  ohannel  aelection  prior  to  regression  retrievals  is  alao 
being  done  in  an  attempt  to  improve  the  above  moisture  results 
particularly.  Only  with  an  increase  in  explained  varianee  will  those 
results  be  acoeptable  as  an  improvement  over  sinple  regression. 


6.0  TRANSMITTANCE  SOFTWARE  FOR  SSB-1  CHANNELS 

Hie  need  for  SSB-1  transmittance  software  arose  ont  of  a  desire  to 
implement  an  iterative  or  minimus- inf creation  retrieval  soheae.  These 
schemes  do  not  require  a  history  of  satellite-RAOB  pairs  which  serve  as 
the  statistical  base  for  regression  or  eigenvector  retrieval  schemes. 
Instead,  they  require  a  knowledge  of  the  atsMspheric  transmittance s  for 
radiation  at  wavelengths  within  the  various  satellite  ehanmel  ranges. 

In  an  iterative  retrieval  scheme  these  tranamittances  can  be  updated 
with  each  adjustment  of  the  temperature  and  moisture  profiles.  In  a 
minimum-information  retrieval  scheme  the  transmittanoes  form  the  basis 
for  the  retrieval  matrix  together  with  a  knowledge  of  the  noise 
variances. 

In  either  retrieval  scheme  tranamittances  need  to  be  calculated 
for  the  appropriate  satellite  channels.  To  meet  this  requirement* 
software  was  adapted  from  various  sources  to  the  S8H-1  channels.  The 
software  employed  polynomial  approximations  of  the  more  precise  line- 
by-line  computations.  Even  with  the  increased  ecaputatiom  speed  which 
results  from  this  approximation*  the  software  still  requires  significant 
computer  resources. 

6.1  UjQ.jjASS^  ChlMfU 

The  SSB-1  HjO  and  filter  response  functions  obtained  from  AF6L 
were  used  to  set  up  a  simplified  polynomial  transmittance 
representation.  The  transmittanoes  calculated  by  this  simplified  and 
fast  method  are  the  basis  for  an  iterative  (physical)  water  vapor  and 
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teaperature  retrieval  scheae. 

Data  froa  three  sources  were  eoabiaed  la  this  endeavor.  Saith 
(1969)  calculated  traasaittaaee  polyaoaials  which  adequately  represented 
the  ataospheric  traasaittaaee  as  it  was  known  at  that  tiae  for  the  15  pa 
C02  and  20  pa  rotational  HjO  bands.  The  Saith  polyaoaials,  as  they  will 
be  called  here,  were  originally  developed  for  the  SHtS  (Satellite 
Infrared  Speetroaeter) .  The  polynoaial  coefficients  are  given  for  every 
5  on-1  froa  about  200  to  820  ea~*.  These  polyaoaials  were  convoluted 
with  the  SSB-1  filter  response  functions  provided  by  AFGL  for  the  HjO 
and  C02  channels.  Then  software  developed  by  NO AA/ NESS,  siailar  to  that 
used  for  window  channel  traasaittaaee  calculations  (Weiareb  and  Hill, 
1980),  was  aodified  to  be  eosqpatiable  with  the  Saith  polyaoaials.  The 
result  was  a  quick  and  easy  tramsaittaaee  calculation  scheae  which  will 
be  explained  in  aore  detail  ae  follows. 

Saith  developed  a  polynoaial  traasaittaaee  representation  as  a 
speedy  alternative  to  the  auch  slower  (but  possibly  aore  aeeurate) 
line-by-line  calculations  which  can  be  used  to  deteraine  the 
transaittance  for  a  given  tsnperature  and  aoisture  profile.  The 
polyaoaials  provided  an  adequate  representation  of  the  eontiauua  and  the 
absorption  bands  as  known  at  that  tiae.  However,  iaproveaents  in 
transaittance  aeasureaents  nay  have  aade  the  polyaoaials  inadequate. 

This  can  be  deteraiaed  by  ecaparison  with  aore  precise  and  up-to-date 
calculations  or  by  eoaparisoa  with  eaperical  transaittance  data. 


Hi*  Sal tli  polynomials  oonsist  of  9  tons  of  the  fon: 
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where  the  Cj's  are  tabulated  transnittance  polynomial  coefficients  at 
every  5  on  1  and  the  A^'s  the  polynomial  tens  given  in  Table  (.1. 
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The  polynomial  representation  was  derived  by  linearising  a 
transmittance  approximation  in  tens  of  optical  mass  of  water  vapor  or 
total  preeipitable  water  (V),  pressure  (P),  and  temperature  (T)  by 
taking  the  logarithm  twice  and  expanding  in  regnlar  polynomials  snoh 
that  in  tens  of  wavenumber  k: 


TAk  *  exp [exp (V) ] 


(6.2) 


The  details  are  given  by  faith  (1969)  or  in  modified  fora  by  Minnie  and 
Cox  (1976). 

In  order  to  specialize  the  polynomials  to  the  sounder  under  study, 
these  polynomials  were  convoluted  with  the  SSB-1  H^O  ad  COj  filter 
response  functions.  Figure  6.1  gives  the  filter  response  for  channel  F3 
on  satellite  VX4537.  The  filter  is  about  30  om”1  wide  and  is 
approximated  every  5  cm-1  (where  the  polynomial  coefficients  are  known) 
by  the  histogram  whioh  gives  an  indication  of  the  quality  of  the 
approximation. 

New  polynomials  were  then  ealeulated  that  were  representations  for 
eaoh  ohannel  by  the  convolution  of  the  tabulated  polynomial  coefficients 
with  the  filter  response  function  fat  eaoh  ohannel. 


where  k  is  the  number  of  faith  polynomials  terms  at  5  cm-1  intervals 
whioh  went  into  the  new  polynomial  term  to  replaee  in  Equation 
6.1.  A  nine-term  polynomial  was.  therefore,  derived  for  each  ohannel  i. 
whioh  would  determine  its  tranamittanoe  for  any  known  pressure, 
temperature,  and  mixing  ratio  profile.  However,  for  the  speoial  case  of 
ohannel  F3  at  536  om  1  and  the  00^  channels,  whioh  experience  absorption 
by  both  EjO  and  COj.  the  separate  ^0  and  COj  tranamittanoe  a  eaoh 
produced  a  polynomial  of  the  form  of  Bquation  6.1.  The  resulting 


FILTER  RESPONSE  (%) 


420  cm*1  Satellite 


Figure  6.1.  Filter  response  function  for  SSH-1  channel  F3.  The  histo¬ 
gram  represents  the  approximation  over  3  cm-1  intervals  for  nse 
with  Smith  transmittance  polynomials. 


transmittance s  then  were  multiplied  together  for  the  total 
transmittance. 


T total  *  *^0  TC02  (6.4) 

The  transmittances  were  then  nnmerically  differenced  in  the 
vertical  to  produce  weighting  functions,  as  shown  in  Figure  6.2  for  the 
H20  channels.  The  weighting  functions  for  all  8  HjO  channels  (F1-F8) 
are  given  along  with  the  separate  H20  and  C02  weighting  funetiona  for 
the  FS  channel.  The  weight  axis  units  are  relative  to  the  pressure 
coordinate  system  used  in  the  vertioal  differencing.  These  weighting 
functions  are  calculated  for  the  1962  US  standard  temperature  profile 
and  for  a  given  moisture  profile  with  a  mixing  ratio  of  lOg  kg"1  at  100 
kPa  (lOOOnb)  and  a  total  preeipitable  water  of  2.55  cm.  These 
temperature  and  moisture  profiles  are  plotted  in  Figures  6.3  and  6.4. 

6.2  CO^  Transmittance  Problem  and  Solution 

The  first  attempt  at  obtaining  the  SSB-1  CO^  channel  transmittances 
proved  to  be  far  from  satisfactory.  Many  of  the  C02  channel  weighting 
functions  were  remarkably  similar  with  maximum  weight  well  below  the 
tropopause  (typically  below  500  mb).  This  was  in  strong  contrast  to  the 
weighting  functions  for  the  Vertical  Temperature  Profile  Radiometer 
( VTPR)  C02  channels  which  have  maxima  spread  throughout  the  stratosphere 
and  troposphere  (MeMillin,  et  al.,  1973).  For  example,  the  weighting 
funetion  for  the  least  transparent  15  pm  channel  (668  cm-1)  should  peak 
well  above  the  tropopause,  aa  was  not  the  case  with  the  polynomial- 
approximated  weighting  funotion. 

These  problems  lead  to  an  evaluation  of  some  software  obtained  in 
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Figure  6.2.  Weighting  functions  for  the  8  SSH-1  BLO  channels  (F1-F8) 

with  separate  weighting  functions  for  the  I,<Taad  CO,  transaittsnee 
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Figure  6.4.  Moisture  profile  with  2.55  ea  of  precipitable  water  weed 
for  traasaittaace  oalealatioae. 


1975  from  NOAA/NESS  for  tho  calculation  of  the  VTPE  channel 
transmittance s.  That  software  was  nsed  for  VTR  data  from  the  NOAA 
satellite  series  and  is  based  on  the  work  of  MoMillin  and  Fleming 
(1976) .  Because  of  the  similarity  between  the  VTPR  and  SSH-1  CO^ 
channels  in  both  central  wavenumber  and  half  width  the  VTHt 
transmittance  software  was  used  as  a  guide  in  interpretng  the  problems 
with  the  software  developed  for  SSH-1  and  eventually  as  replacement 
software  for  the  00^  channel  transmittanoes.  A  comparison  of  the  VTTR 
channels  with  the  corresponding  SSH-1  channel s  is  given  in  Table  6.2. 

Table  6.2 

Comparison  of  VCT  and  Equivalent  SSH-1  Channels 
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4 
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10. 

E5 

725. 

10. 
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10. 
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535. 

18. 

F5(E7) 

535. 

16.5 

Window 

8 

833. 

10. 

W2(E8) 

835. 
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The  problem  with  the  transmittanoes  for  the  SSH-1  (X>2  channels 
(E1-E6)  seemed  to  stem  from  an  inoorrect  application  of  the  Smith 
polynomials  to  the  CO^  absorption  in  the  15  pm  region.  To  remedy  this, 
the  C02  absorption  part  of  the  VTPE.  transmittance s  was  adopted. 

In  comparison  of  the  two  sets  of  transmittance  software  another 
possible  error  was  discovered.  The  discrepancy  was  a  variation  in  the 
definition  of  the  'X'  term  of  Table  6.1.  By  a  comparison  of  the 
polynomials  generated  here  for  the  SSH-1  ohannels  and  the  polynomials 


given  for  the  ViVk  H^O  absorption,  a  modification  was  made  to  the 
software  to  make  the  calculations  consistent.  Surprisingly  the  results 
were  not  affected  very  severly  by  this  change. 

The  final  CO^  weighting  functions  are  given  in  Figure  6.S.  Again, 
these  transmittance s  are  specific  to  the  temperature  and  moisture 
profiles  used  in  the  calculations.  Before  extensive  use,  the 
transmittance  software  should  be  compared  to  line-by-line  calculations 
or  to  experimentally  observed  transmittences.  If  the  transmittance s  are 
not  good  representations,  bias  errors  in  retrievals  could  result,  but  at 
least  the  transmittanoes  should  be  known  well  enough  in  order  to 
adequately  reflect  the  temperature  and  moisture  changes  necessary  in  the 
retrieval  process.  This  polynomial  transmittance  representation  can 
then  be  used  to  calculate  the  transmittance s  and  resulting  weighting 
functions  which  are  needed  in  a  physical  retrieval  scheme. 

6.3  WAn.4.o.iL_<?h.W9l 

Finally,  transmittance  software  based  on  the  work  of  Weinreb  and 
Hill  (1980)  was  used  for  the  12  pm  window  channel  on  the  SSH-1 
instruments.  This  software  was  originally  developed  for  broadband 


window  channels  on  the  TIROS-N/NOAA  satellite  series 
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Figure  6.5.  Weighting  functions  for  SSH-1  CO.  channels  (E1-E6) 
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7 .0  SUMMARY 


Several  aspects  of  the  infrared  remote  sounding  problem  particular 
to  SSH-1  have  been  addressed.  Each  of  theae  areas  of  research  is 
related  to  the  two  objectives  listed  in  the  introduction.  For  research 
purposes  real  (not  simulated)  SSB-1  radiance  data  with  corresponding 
rawinsonde  (RAOB)  data  were  obtained  from  AFGVC.  These  SSH-l/RAOB  pairs 
covered  the  period  from  June  1979  through  February  1980. 

A  structure  analysis  was  performed  on  the  SSH-1  radiances  in  order 
to  estisiate  the  radiance  noise  levels.  The  results  confirm  that  the 
SSH-1  radiances  have  not  been  degraded  significantly,  if  at  all.  from 
their  expected  noise  levels.  (These  noise  levels  may  be  translated  into 
the  uncertainty  to  be  expected  for  temperature  and  moisture  parameters 
provided  by  perfeot  retrieval  schemes  using  this  noisy  data.)  In  other 
words,  the  problems  with  SSH-1  retrievals  are  not  caused  by  noisy 
radiances.  Eivenveetor  analysis,  on  the  other  hand,  does  indicate 
significant  redundancy  in  the  ohannels  which  limits  the  amount  of 
information  available. 

The  scan-angle  corrections  as  applied  by  AFGVC  were  found  to  be 
adequate  in  all  but  one  channel  (F8).  The  solution  was  to  apply  the 
correction  polynomial  for  a  similar  ehannel  (FI)  to  the  one  channel  (F8) 
for  which  the  polynomial  produced  a  large  bias.  These  scan  angle 
corrections  are  necessary  for  the  application  of  a  statistical  retrieval 
scheme  which  requires  radiances  from  a  eommon  nadir  (zenith-angle 
corrected)  viewidg  position. 


The  cloud  detection  problen  ie  also  very  important  ainoe  this  is 
the  source  of  sumy  problems  with  infrared  retrieval  seheaes.  Snail 
amounts  of  clond  as  veil  as  cirms  are  very  hard  to  detect  and 
eliainate.  Two  nethods  are  presented;  one  which  depends  on  radiance 
thresholds  which  are  deterained  enperically  or  theoretically;  and  the 
other  nethod  which  can  be  applied  to  future  SSH-2  data.  The  second 
nethod  relies  on  the  use  of  two  window  channels  with  different  spectral 
characteristics.  Both  cloud  aaount  and  teaperature  (or  height)  can  be 
deterained  for  a  single  FOV  if  a  surface  or  background  teaperature  is 
known  or  assuaed.  This  nethod,  however,  assuaes  that  the  field-of~view 
is  only  partly  cloudy  and  contains  a  unifora  eloud  top.  Even  if  these 
assumptions  are  not  net,  the  two  nethods  can  effectively  detect  cloudy 
fields-of-view  which  can  then  be  eliminated  from  consideration. 

The  assessaent  of  a  regression  retrieval  scheae  for  teaperatures 
and  aoisture  does  indicate  that  the  SSH-1  radiances  do  contain 
teaperature  and  aoisture  information.  Results  of  retrievals  on 
independent  data  (data  other  than  that  used  to  develop  the  regression 
coefficients)  are  degraded  from  dependent  data  tests  but  are  not 
ooapletely  noise.  Moisture,  especially  of  the  lowest  layers,  oan  be 
retrieved,  contrary  to  prior  experience  with  SSH-1  data  at  AFGWC.  These 
results  are  also  oonfiraed  by  others  who  have  worked  with  siaulated  data 
for  SSB-1.  Deficiencies  in  the  statistical  data  base  used  by  AFGWC  for 
the  regression  analysis  aay  be  the  source  of  the  problem. 

Finally,  transaittanee  software  for  the  SSH-1  H^O,  CO^,  and  window 
channels  was  developed.  The  need  for  this  software  arose  froa  the 
expected  application  of  an  iterative  aoisture  and  teaperature  retrieval 
algoritha.  However,  because  of  computing  expense  and  tine  liaitations 


tills  course  of  research  was  discontinued.  Much  swifter  regression 
retrievals  will  most  likely  he  adopted  with  no  degradation  of  retrieval 
quality.  Appropriate  selection  of  the  SSH-2  channels  whioh  are  to  be 
used  in  retrieval  will  be  neoessary. 

Much  more  work  is  needed,  especially  to  obtain  useful  moiature 
information.  The  moisture  channel  radianoea  seem  to  be  performing  as 
expected,  but  the  channel  selection  is  far  from  optimised  since  the 
eigen-analysis  shows  that  there  are  fewer  than  half  as  many  independent 
pieces  of  information  than  there  are  SSB-1  radiance  channels.  However, 
some  combination  of  SSH  channels  should  produoe  limited  moisture 
information,  especially  when  need  with  the  two  window  channela  on  the 
upcoming  SSH-2.  The  most  practically  retrievable  quantity  may  prove  to 
be  integrated  moisture,  with  little  or  no  vertical  discrimination. 
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